I  REPORT  DOCUMENTATION  PAGE 

i  ’o«f  ow  •««»•» 

Im  A  f  ^  j  I  3*«t»n9  400  '«.•#<*. i*g  :»*  :H«ct>on  of  m»o^*"«li' 

*  *  *  ®  *■  "  m  "Qucirq  th<t  Ourd*''  'C  A*Sh>r»atOT>  H*i(jaujrt#f 


form  Approve 

OMB  Vo  0704  0198 


TTT 

III 

r 

i  ii! 

F 

pi  r 

jlFf] 

ill 

li 

iJ 

L 

'll  1 

Wii 

m 

l>C*  '  *tr.«4r#0  *J  •’uu'  o*y  MOOrw  >n<*udirg  tne  tm«  for  r*v>9m<nq  ft*i^rt  t'-s  ,««rc*.»3  »«TTr^T!^*^^r^ 

3»«tang  «*.  #**,i*<5  :»*  cittcron  of  mfoT'Y'Btio*  Send  commwti  rM*ro«*q  thi*  t>ur<3«o  •M.F**re  ,t  «r-,  .»»-*?  j  w,  ... 
■duonq  tht»  ouffle'*  •■?  A*$h.p»qion  Hf*aou*rref%  ServcM.  Director*!*  ?o r  ■  nfop»w*t,on  Ooe'*;-on\  *~<j  a^K(.s  t  ^  .„ 

And  tc  !-•  O**  ■*  •  M.*n4q*m**t  ap*  Budget  P-ipef won*  **duCTionproi«<UO >0*  0 ’00)  A*sn  nqicn  C2CSCJ 


2  REPORT  OATC 

April  22,  1994 


S.  FUNDING  NUMBERS 


3.  REPORT  TYPE  AND  OATES  COVEREO 

TECHNICAL  REPORT 


Isolated  S042-  and  PO43-  Anions  Do  Not  Kxist 


6  AUTHOR(S) 

Alexander  I.  Boldryev  and  Jack  Simons 


G  N00014-89-J-L497 
R&T  4131050 


7  PERFORMING  ORGANIZATION  NAME(S)  AND  AOORESS(ES)  i 

UNIVERSITY  OF  UTAH 
DEPARTMENT  OF  CHEMISTRY 
SALT  LAKE  CITY,  UTAH  84112 
U.S.A. 


kS 


D  i  iy 

ELECTE  I 

»KX04t9W 


9.  SPONSORING  .  MONITORING  AGENCY  NAME(S)  AND  ADDRE 

OFFICE  OF  NAVAL  RESEARCH 
CHEMISTRY  PROGRAM 
800  NORTH  QUINCY  ST. 

ARLINGTON,  VIRGINIA  22217-5000 


11.  SUPPLEMENTARY  NOTES 

J.  Phys.  Chem.,  98,  2298  (1994) 


12«.  DISTRIBUTION  AVAILABILITY  STATEMENT 


THIS  DOCUMENT  HAS  BEEN  APPROVED  FOR  PUBLIC  RELEASE  AND 
SALE;  ITS  DISTRIBUTION  IS  UNLIMITED 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


Technical  Report  No.  42 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 


12b.  DISTRIBUTION  CODE 

UNLIMITED 


13.  ABSTRACT  (Maximum  200  word*) 

Ab  initio  theoretical  methods  that  incorporate  effects  of  orbital  relaxation  and  electron  correlation  are  used 
to  examine  the  possibility  that  doubly  or  triply  charged  oxy-anions  of  the  class  SO4-2,  SeO^2  TeC>4'  2  or 
P04"3.  ASO4'3,  SbC>4‘3  are  stable  relative  to  electron  detachment  in  the  absence  of  any  surrounding  medium. 
In  all  of  these  cases,  our  findings  indicate  quite  clearly  that  these  multiply  charged  oxyanions  are  not  stable. 


15.  NUMBER  OF  PAGES 


16.  PRICE  COOE 


17.  SECURITY  CLASSIFICATION  118,  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIFICATION  20.  LIMITATION  OF  ABSTRACT 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  '  UNCLASSIFIED  UNCLASSIFIED 


UNCLASSIFIED 


Standard  form  298  2  89) 

*'WCT.to*d  bv  amsi  Std  j)9  9 
290-102 


I 


OFFICE  OF  NAVAL  RESEARCH 

Contract  N00014-89-J-1497 

R&T  Code  4131050 
Scientific  Officer:  Dr.  R.  Nowak 

Technical  Report  No.  42 

"Isolated  S042'  and  P043*  Anions  Do  Not  Exist” 
,  by 

Alexander  I.  Boldyrev  and  Jack  Simons 
Prepared  for  Publication  in 
The  Journal  of  Physical  Chemistry 

The  University  of  Utah 
Department  of  Chemistry 
Salt  Lake  City,  Utah  84112-1194 

April  22,  1994 


Reproduction  in  whole  or  in  part  is  permitted  for  any 
purpose  of  the  United  States  Government 

This  document  has  been  approved  for  public  release  and 
sale;  its  distribution  is  unlimited. 


f  (  H-  t  (') 
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Abstract 


Ab  initio  theoretical  methods  that  incorporate  effects  of  orbital  relaxation 
and  electron  correlation  are  used  to  examine  the  possibility  that  doubly  or  triply 
charged  oxy-anions  of  the  class  SCV2,  SeO 4*2,  TeCV2,  or  PO4-3,  ASO4'3, 
Sb04’3  are  stable  relative  to  electron  detachment  in  the  absence  of  any 
surrounding  medium.  In  all  of  these  cases,  our  findings  indicate  quite  clearly 
that  these  multiply  charged  oxyanions  are  not  stable. 
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*  I.  Introduction 

Multiply  charged  anions  such  as  SO4'2  and  PO4-3  play  important  roles  in 
chemistry  and  biochemistry,  where  they  appear  as  structural  units  of  crystals, 
solutions  and  molten  substances.  The  bond  lengths1  and  vibrational 
frequencies2,  effective  ionic  radii  and  electronic  polarizabilities3  as  well  as 
thermodynamic  functions4  of  SO4'2  and  the  bond  lengths  of  PCV3-5  have  been 
determined  experimentally.  However,  all  of  these  data  were  obtained  in  some 
medium,  which  may  act  to  stabilize  and  hence  change  the  properties  of  the  free 
anions.  The  question,  "Can  these  free  anions  exist?”  is  therefore  still  open. 

The  experimental  results  for  free  SO4*2  are  controversial. 
Cancerevicius63  and  Schulten  and  Rollgen6b  reported  observation  free  SO4'2 
using  field  desorption  mass  spectrometry  of  inorganic  ions  from  a  polymeric 
matrix.  However,  several  others  experimental  studies7  were  not  able  to  detect 
SO4-2  in  the  gas  phase.  In  theoretical  work,  Scheller  and  Cederbaum8d  found 
that  CO3-  and  a  free  electron  are  more  stable  than  CO3-2,  so  free  CO3-2  does 
not  exist  Recently,  doubly  charged  anions  such  as  BeF4-2,8a  SeFs'2,  TeCIs'2, 
TeF8'2,8b  LX3-2  (L=Li,Na,K;  XsF.CI)80-11  have  been  studied  using  ab  initio 
theoretical  methods  where  it  has  been  shown  that,  although  all  of  these  species 
are  not  stable  thermodynamically  with  respect  to  loss  of  one  F*.  they  are 
geometrically  metastable  and  electronically  stable  species  with  large 
dissociation  barriers  and  hence  large  lifetimes. 

II.  Methods  Used 

In  this  article  we  present  results  of  our  calculations  for  SO4'2  and  PO4'3  , 
as  well  as  for  the  valence  isoelectronic  Se04'2,  TeC>4'2,  ASO4-3  and  SbC>4'3 
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ions.  The  bond  lengths  of  SO4'2  and  P(V3  were  optimized  by  employing 
analytical  gradients9  using  a  polarized  split-valence  basis  set  (6-31 1+G*)10  at 
the  MP2(full)  level.  The  MP2  (full)/6-31+G*  equilibrium  geometries  were  then 
used  to  evaluate  electron  correlation  corrections  in  the  frozen-core 
approximation  by  Moller-Plesset  perturbation  theory  to  full  fourth  order11  and  by 
the  (U)QCISD(T)  method12  using  6-31 1+G*  basis  sets.  The  UHF  wave  functions 
for  open  shell  systems  were  projected  to  pure  spectroscopic  states  and  the 
corresponding  results  are  denoted  PUHF,  PMP2,  PMP3  and  PMP4.13  Bond 
lengths  for  the  isoelectronic  SeCV2,  TeCV2,  AsCV3  and  Sb04*3  have  been 
obtained  using  the  MP2  approximation  and  the  Los  Alamos 
pseudopotentials.14  We  used  Dunning's  valence  double-zeta  basis  set 
(9s5p/3s2p)  for  the  oxygen  atoms15  and  valence  2s2p  basis  sets  for  the 
selenium,  tellurium,  arsenic  and  antimony  atoms14  extended  by  polarization  d- 
functions  on  all  atoms  and  diffuse  s-  and  p-functions  on  the  oxygen  atoms. 
Exponents  for  the  polarization  d-functions  (0.293  (As),  0.338  (Se),  0.211  (Sb), 
0.237  (T e))  were  taken  from  Ref.  16  and  exponents  of  the  diffuse  s-  and  p- 
functions  0.059  were  taken  from  Ref.  17.  Finally  we  calculated  the  total  energies 
of  SO4'2,  Se04'2,  TeCV2  and  the  corresponding  singly  charged  anions  (with 
the  detachment  of  one  electron  from  the  Iti-HOMO  and  from  the  next  lowest 
3t2*HOMO-l)  at  the  optimal  geometry  of  the  doubly  charged  anions  at  the  MPn 
(n=2-4),  QCISD,  QCISD(T)  levels  and  using  the  recently  suggested  Feenberg 
perturbation  series  (For  details  see  Ref.  18  and  references  therein)  FPTn  (n=2- 
4).  The  Gaussian  92  program  suite19  was  used  to  perform  all  of  the  calculations 
whose  results  are  discussed  here. 
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III.  Results 

Our  optimal  S-0  bond  length  1 .52  A  for  isolated  SO42  is  close  to  the 
corresponding  values  1.47-1.50  A  in  various  crystals,  1.50  A  in  molten  salts  and 
1.50  A  in  solutions.  The  orbital  energy  of  the  Iti-HOMO  is  negative  (e=-0.73  eV) 
which  means  that  this  dianion  is  electronically  stable  within  the  Koopmans' 
theorem  approximation.  Because  the  Koopmans'  level  description  involves 
doubly  occupied  molecular  orbitals,  it  is  unable  to  describe  electron  ejection  to 
SO4-  +  e‘.  By  taking  orbital  relaxation  into  account  (i.e.,  carrying  out  separate 
SCF  level  calculations  on  SO4*2  and  SO4*)  in  the  so-called  ASCF  treatment,  the 
total  energy  of  SO4'  is  found  to  be  1.10  eV  lower  than  SO4-2  at  the  optimal 
dianion  geometry.  Hence,  the  doubly  charged  anion  is  predicted  to  be 
electronically  unstable  at  the  ASCF  level.  Because  differential  electron 
correlation  energies  can  amount  to  leV  or  more,  it  is  therefore  essential  to 
examine  the  dianion  and  anion  stabilities  at  various  correlated  levels.  Doing 
so,  we  find  that  correlation  corrections  to  the  ASCF  prediction  calculated  at  the 
PMP2,  PMP3,  PMP4  and  QCISD(T)  levels  with  6-31 1+G*  basis  sets  increase 
the  instability  of  SO4'2  with  respect  to  SO4'  by  approximately  0.2-0.4  eV.  At  our 
highest  PMP4  and  QCISD(T)  levels,  the  dianion  is  less  stable  than  the  singly 
charged  anion  by  1.57  eV  and  1 .34  eV,  respectively  (see  Table  1).  Because  the 
final  state  of  the  singly  charged  SCV  anion  in  the  2Ti  electronic  state  is  not 
geometrically  stable  at  T<j  symmetry  due  to  Jahn-Teller  distortion,  the  SO4' 
anion  will  be  even  more  stable  at  its  optimal  geometry,  therefore  SO4'2  is  even 
more  unstable  with  respect  to  electron  detachment.  Weikert  and  Cederbaum8® 
also  have  found  that  S042*  dianion  is  unstable  to  electron  loss  according  to  the 
Green's  function  calculations. 
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The  next  occupied  3t2-(HOMO-1)  orbital  is  more  stable  than  Iti-HOMO 
by  1 .8  eV  at  Koopmans'  approximation,  however  due  to  a  larger  electron 
relaxation  effect,  at  the  ASCF  level,  the  energy  difference  between  these  two 
states  is  only  0.7  eV  with  the  electronic  state  lying  lower.  When  electron 
correlation  energy  is  taken  into  account,  the  ^T^  state  of  SCV  is  only  0.2-0. 4  eV 
lower  than  the  2T2  state  at  the  optimal  geometry  for  SCV2.  Calculated  vertical 
detachment  energies  for  both  states  at  the  FPT4  level  are  in  the  range  of  the 
corresponding  values  at  the  PMP4  and  QCISD(T)  levels.  Because,  at  all 
correlated  levels,  SO4'2  is  not  electronically  stable,  our  final  conclusion  is  that 
free  SO4'2  does  not  exist  as  a  stable  species,  which  agrees  with  most  of  the 
available  experimental  data.7  Our  findings  do  not  preclude  the  possibility  that 
SO4'2  can  exist  as  a  short-lived  resonance  state,  but  there  exists  in  our  data  no 
hint  of  such  behavior. 

Based  on  previous  calculations,20  it  is  known  that  when  the  distances 
between  the  electronegative  ligands  in  closed  shell  multiply  charged  anions 
increase,  the  stability  of  the  HOMO  increases  also.  Therefore,  we  expect  the 
HOMOs  of  both  Se04  2  and  Te04  2  to  be  more  stable  than  that  of  SO42  and 
perhaps  one  of  these  dianions  can  survive  in  the  free  state.  Indeed,  the  orbital 
energies  of  the  HOMOs  varies  in  this  manner:  SO4'2  (-0.73  eV),  Se04'2  (-1  55 
eV),  Te04*2  (-2.1 1  eV)  and  at  the  Koopman's  level  both  dianions  are  very  stable 
with  respect  to  electron  loss. 

The  optimal  bond  length  R(Se-0)=1 .678  A  of  free  Se04*2  is  very  close  to 
the  corresponding  values  in  crystal  salts:  1.66±0.01  A  ((NH4)2S0421a), 
1.654±0.021  A  (Na2Se0421b),  1.65±0.01  A  (K2Se0421c)  and  1.64  A 
(Rb2Se0421d).  The  corresponding  Te  crystal  salts  Li2Te0422a  and  Na2Te0422b 
are  not  composed  of  Te04'2  anions.  Instead,  the  tellurium  atoms  occur  in 
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octahedral  oxygen  coordination  in  these  salts.  Our  optimal  bond  length  for 
Te04'2  R(Te-O)  is  equal  to  1 .839  A. 

Both  Se04'2  and  Te04'2  have  local  minima  at  the  tetrahedral  structure 
because  all  of  the  computed  vibrational  frequencies  are  positive(Se04'2: 
vi(ai)=802  cm*1,  V2(e)=310  cm-1,  V3(t2)=819  cm'1,  V4(t2)=405  cm*1  and  Te042: 
vi(ai)=738  cm'1,  V2(e)=250  cm'1,  V3(t2)=744  cm*1,  V4(t2)=318  cm'1,  all  at  MP2). 
Therefore,  if  we  keep  the  number  of  electrons  constant,  these  dianions  are 
geometrically  stable.  However,  when  electron  relaxation  and  electron 
correlation  are  taken  into  account,  both  Se04'2  and  Te04'2  are  not 
electronically  stable  with  respect  to  removal  of  an  electron  from  the  Iti-HOMO 
to  produce  a  2Ti  final  electronic  state  (see  Table  1).  Therefore,  we  conclude 
that  both  of  these  free  doubly  charged  anions  do  not  exist  either. 

An  important  consequence  of  these  calculations  is  that  even  if  a  multiply 
charged  anion  is  very  stable  at  the  Koopmans1  approximation  level  (as  happens 
for  TeCV2),  we  should  not  trust  these  results.  Electron  relaxation  and  electron 
correlation  are  extremely  important  even  for  making  qualitative  predictions. 

Based  on  the  results  presented  above  for  doubly  charged  anions,  we  did 
not  expect  that  any  of  the  triply  charged  PO4'3,  AsCV3  and  SbCV3  ions  could 
survive  even  as  metastable  species.  Our  computational  results  confirm  these 
predictions;  all  three  anions  have  very  positive  HOMO  energies:  +5.33  eV  (PO4' 
3),  +4.33  eV  (ASO4'3)  and  +3.52  eV  (Sb04'3)  at  the  corresponding  optimal 
Td(1  Ai )  geometries  with  bond  lengths  R(P-0}=1.599  A,  R(As-0)=1 .750  A  and 
R(Sb-0)=1 .936  A,  calculated  at  the  MP2/LANL1DZ+  dTe,o+diffs,po  level.  Not 
only  is  the  HOMO  energy  positive,  but  the  energies  of  a  few  other  occupied 
orbitals  are  also  positive.  Therefore,  these  free  triply  charged  anions  are  not 
electronically  stable.  Our  findings  do  not  preclude  the  possibility  that  SO4-2, 
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Se04'2.  T©04-2  and  triply  PO4'3,  ASO4'3,  Sb04‘3can  exist  as  a  short-lived 
resonance  state,  but  there  exists  in  our  data  no  hint  of  such  behavior. 

The  instability  of  the  free  doubly  SO4'2,  Se04-2,  TeCV2  and  triply  P04-3, 
ASO4'3,  Sb04*3  charged  studied  here  has  important  consequences.  For 
example,  in  the  thermochemic  calculations  of  crystal  energies  for  salts 
composed  from  these  anions,  one  should  not  consider  the  anions  as  stable 
structural  fragments,  the  stable  free  fragment  is  actually  the  singly  charged 
anion.  Moreover,  the  structural  parameters,  harmonic  frequencies,  force 
constants  and  other  molecular  properties  for  these  doubly  and  triply  charged 
anions  should  not  be  attributed  to  the  free  anions,  but  to  anions  in  the  a  crystal 
field  that  acts  to  stability  the  electronic  state  of  the  ion. 
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